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ABSTRACT 

One-dimensional, electrostatic, particle-in-cell simulations are used to explore two 
mechanisms proposed to explain turbulent broadening of soft X-ray emission lines o heavy 
ions observed during solar flares and the presence of blue-shifted components. Resu 
the simulations are in qualitative agreement with the observations. 

1. INTRODUCTION , 

Observations of soft X-ray emission lines of heavy ions during solar flares show the 

presence of energetic turbulent motions in the heavy ions and blue-shifted componen s 
(Antonucci et al. 1982; MacNeice et al. 1982; Tanaka 1987). The source of the energy for 
these phenomena is believed to be the energetic beam electrons produced m the impulsive 
phase of a flare. The exact mechanism (or mechanisms) that produces the emission ^ line 
features is a subject of debate. Two such proposed mechanisms, similar m that both utilize 
superheated electron populations to drive the ion heating, are examined using electrostatic 
particle-in-cell simulations. . , , , 

The first mechanism is in-situ coronal maser heating (Fig. ), in w 1C e e 
cyclotron maser radiation emitted by the energetic electrons is absorbed by ambient coronal 
electrons (Melrose and Dulk 1982, 1984; McKean, Winglee and Dulk 1989). The heated 
electrons expand up magnetic field lines in the presence of local density gradients, creat ng 
an ambipolar electric field due to the charge separartion between the electrons and the 
source region ions which then accelerates the ions. It is found that the ambipolau Md 
produces differential ion motion between ion species; this differential motion is unstable 
to the ion-ion streaming instability, which then preferentially heats and accelerates the 
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Fig 1 : Schematic of coronal maser heating mechanism. 
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Fig. 2: Schematic of chromospheric heating mechanism. 
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heavier ion species. The resulting ion distribution functions bear similarities to observed 
soft X-ray line profiles. 

The second mechanism is chromospheric heating (Fig. 2), in which chromospheric 
electrons are collisionally heated by energetic electrons (Winglee 1989). Under certain con- 
ditions, the collisional heating of the chromospheric electrons exceeds the collisional heating 
of the ions (Winglee 1989); the electrons then expand up field lines into the corona, which 
results in the formation of an ambipolar electric field and resultant ion acceleration. Heat- 
ing due to the ion-ion instability is less important in this case, though turbulent-appearing 
spectra may still be produced due to spatially inhomogeneous ion velocity profiles. 

2. SIMULATION MODEL 

One-dimensional, electrostatic, particle-in-cell simulations are used in this paper to 
model the self-consistent evolution of plasma during the outward expansion of a super- 
heated electron population into an ambient plasma at higher altitude and lower density, 
including the development of an ambipolar electric field and of wave-particle interactions. 

The plasma is assumed to have two ion species. One species of ion with mass mjj 
corresponds to H + , the dominant species in the solar atmosphere, and accounts for 90% of 
the total ion number in these simulations. The second ion has a mass mj e = 4 m#, so it 
is similar in behavior to He + or any other ion with q/m — 1/4. These heavier ions have a 
charge-to-mass ratio that roughly corresponds to an ion that has yet to be fully stripped of 
its electrons. The reason for simulating an ion species that has not yet been fully stripped 
is that, at the time the ions are being accelerated, they have not been in contact with a 
population of hot, energetic electrons long enough to have reached an equilibrium, fully 
stripped ionized state. In order to allow ions to react to electron motion on the time scale 
of the simulation, it is assumed that mu = 100m e , where m e is the electron mass. 

The length of the simulation system (Figure 3) is 2048A, where A is the simulation 
cell size. The system is divided into three sections. Between x = 0 and x/A — 384 lies 
Region A, which is a region of enhanced plasma density «*• This region represents the 
lower corona in the coronal heating simulation and the chromosphere in the chromospheric 
heating simulation. Between x/A = 512 and x/A = 2048 lies Region C, which represents 
the ambient coronal plasma where the plasma density nc is less than in Region A. Between 
x/A = 384 and x/A — 512 lies Region B, where there is a linear change in density between 
Regions A and C. The enhanced density in Regions A and B, with a hot electron population 
that has a temperature elevated significantly above the ion temperature, is the ultimate 
source of the free energy available for the formation of the ambipolar electric field and 
for the growth of the various instabilities seen in these simulations. Therefore, particles 
that begin the simulation in Regions A and B are referred to collectively “source region 
particles.” In the chromospheric heating simulation, there is also an ion temperature 
gradient across Region B, with the ion temperature higher in the corona than in the 
chromosphere. 
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The plasma in these simulations is assumed to be collisionless. Particles are reflected 
after crossing either boundary of the system. 


Region R Region B Region C 



X /A 


Fig. 3: Schematic model of simulation system showing where 
Regions A, B, and C lie and the initial ion and electron 
temperature and density profiles in the two simulations 
whose results are shown in this paper. 


3. IN- SITU CORONAL HEATING 

The effects of extended in- situ coronal heating are examined in this section. The 
electrons are assumed to be uniformly heated with the temperature uniform through the 
system and elevated above the ion temperature. Due to the assumed density gradient, there 
is a pressure excess in the dense region which drives an outward propagating conduction 
front and subsequent plasma heating. The plasma density changes by a factor of four 
across Region B, i.e., n^/nc = 4. The ratio of electron temperature to ion temperature 
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T t /T{ = 100. The acoustic velocity c 9 — (Te/m#) 1 / 2 = 1.18 a/# A, where uh is the H+ 
plasma frequency in Region A. 

The electrons expanding out of the source region set up an ambipolar electric field 
due to the charge separation between the electrons and the source region ions. The po- 
tential drop across the ambipolar field can be inferred from the maxmimum velocity of 
the reflected electrons, which is about 2 vr«> implying a potential drop of ~ 4 T e . This 
potential drop is strong enough to confine the bulk (~ 95%) of the source region electrons 
to the region within and behind the ambipolar field. 

The ion dynamics is illustrated in Figure 4, which shows the development of the ion 
phase space, with the H+ ions in the left column and the He+ ions in the right column. 
The initial phase spaces are shown in Figs. 4a and 4b. Once formed, the ambipolar field 
accelerates both ion species outwards in bulk (Figs. 4c and d). The acceleration of the if + 
ions to high velocities occurs before the He+ ions because of their higher charge-to-mass 
ratio; thus the average velocity of the accelerated if + is ~ 4 times higher than the same 
for the accelerated ife+ at t out = 30. 

The ambipolar field, initially limited in spatial extent to Region B, expands beyond 
that region as the footpoints of the field (those points at which the H + bulk velocity equals 
zero) propagate away from Region B at about the local acoustic velocity c 9 . By u;//£ — 30, 
the ambipolar field extends from x/A ~ 340 to x/A ~ 560 (Fig. 4c). 

The footpoints represent the movement of the two fronts of the ambipolar field, 
one of which (a compression front) propagates into and compresses the ambient ions, 
accelerating them forward; the other (an expansion front) moves into and expands the 
source region ions, drawing them outward. The two fronts initially overlap; when they 
separate, a velocity plateau forms in the H+ phase space (Fig. 4e), and the electric field 
is zero or near zero in the region of the plateau. 

After the fronts separate, the expansion front (marked by the right side dashed line 
in figures 4c-l) propagates to the left at ~ 0.65 c 9 . It constantly expands in size with time 
because its footpoint moves inward at a higher velocity (~ c*) than the entire front does. 
By contrast, the size of the compression front (which is marked by the left side dashed 
line in Figures 4c-l) front is monotonically reduced (Fig. 4g) until it becomes a narrow, 
shock-like structure (Fig. 4i, k), with the footpoint of the front moving to the left at ~ c 9 
and the entire front at ~ 1.3 c 9 . The front does not break because, as it steepens, the 
dispersion of its energy increases and effectively counters the steepening. 

By the times in in Figures 4i and k, the front can be described a s a shock because 
there are sharp discontinuities in bulk plasma properties; in other words, large changes in 
plasma properties such as density and bulk velocity take place across a very narrow region 
(on the order of an ion acoustic wave wavelength wide). This shock is important because 
it marks the leading edge of a density wave associated with the expansion of the source 
region plasma. Its movement through a coronal loop may be seen in spatially resolved 
X-ray observations because the shock should leave an increased soft X-ray emission in its 
wake due to the sharply increased numbers of energetic electrons that are confined behind 
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the shock. Moreover, since it is travelling only at ~ c*, high time resolution is not required 
to resolve its propagation. 

The compression front bears similarities to the ion sound conduction front of Brown, 
Melrose, and Spicer (1979) in that it confines energetic source region electrons to the region 
behind the front and reduces the heat flux outwards to a value determined by the velocity 
of the front, which is of order ~ c e . This limitation of the heat flux is the principle 
function of a conduction front, which means that the compression front can be identified 
as a conduction front as well. Because of the different confinement mechanism, it is best 
described as an “ambipolar field conduction front” to differentiate it from the ion sound 
conduction front of Brown, Melrose, and Spicer (1979). 

The passage of the two fronts through the plasma accelerates H+ ions to a bulk 
velocity, hereafter the “beam” velocity, ~ 0.65 c,, which corresponds to a “beam” energy 
of ~ 0.42 T e . The 2J + ions clearly do not experience the potential drop of ~ 4 T t that the 
electrons do. The explanation for this apparent discrepancy is that the ambipolar field 
fronts move through the ions in less time than it would take for an ion to experience the 
full potential drop if the fronts were stationary. Because of their smaller charge-to-mass 
ratio and hence mobility, He* ions gain less energy than H+ ions in the potential drop of 
a front. 

Ion acoustic waves excited by expanding heated and return current electrons are 
manifested as high energy tails in the H + ion phase space (Fig. 4e) but not the He+ phase 
space (Fig. 4f). These tail ions can reach energies as high as T e . The highest energy H+ 
ions are produced in and close to the conduction front (Figs. 4e and g). This phenomena is 
a result both of the ambipolar field and of the rapid dissipation of ion waves in this region, 
which happens because the steep velocity gradient increases the effective temperature of 
the H+ ions in the front, thereby reducing T e /T{. The dissipated wave energy goes into 
the production of a population of high energy ions, which are further accelerated through 
interaction with the ambipolar field. These ions, hereafter called “fast beam ions,” can 
reach velocities as high as ~ 2.6 c SJ or energies of ~ 7 T e = 700 T, (Figs. 4g, i and k). 

The differential motions between the H + ions and the He+ ions cam excite an ion- 
ion streaming instability which transfers energy from the faster to the slower ions. The 
requirements for this instability to appear in a one-dimensional plasma are (a) that the 
ratio T t /Ti > 1 and (b) that c s > A u, > where Au, = u H - u H t is the difference 
between the H+ and Ht + bulk velocities and v* < c $ is the threshhold velocity for the 
instability (Dusenberry et al. 1988; Winglee 1989). Two important signatures of the 
ion-ion instability axe: (i) nonthermal heating and acceleration of the heavy ions, which 
produces a distribution function with a broad plateau between the initial H+ and Hc+ 
bulk velocities and a corresponding increase in the He + bulk velocity (Winglee et al. 1989); 
and (ii) the formation of a large, nonthermal low-energy tail on the H + distribution and 
a corresponding decrease in the H + bulk velocity (Winglee et al. 1989; Winglee 1989). 

These signatures of the ion-ion instability can be seen at later times in Figure 4. 
Heavy ions in the beam region and in the fronts have undergone some heating by = 150 
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(Fig. 4h); the velocity dispersion increases further over the remainder of the run (Figs. 4j 
and 1). The heating causes an increase in the maximum velocities of the He + ions to the 
point where they are comparable to the H+ beam velocity. In the phase space plots, 
a low energy tail develops between v = 0 and v = t>t in the beam region (Fig. 4g). This 
low energy tail increases in density with time (Figs. 4i, k). 

Figures 5a and b shows the variation of the density with x at = 0, 150 and 
300 for the H + and He + ions. The bulk flow of H + ions outward reduces the density in 
the source region, though the simulation ends before any significant depletion occurs at 
the left boundary. The conduction front is characterized by a sharp jump in H + density, 
which reflects the shock formation at that front. Behind the front, the H + ions have a 
uniform density between x/A — 450 and 650 at wj/t = 150, and x/A — 450 and 900 at 
u/j jt = 300. Thus, in this region, the ambipolar field has destroyed the density gradient, 
and, since this region corresponds to the velocity plateau in Figures 4g and k, its fronts 
have left the region. 

The He + density profile (Fig. 5b) lags the H+ profile because of the slower bulk 
movement of these ions. This lag results in an enhancement of He+ with respect to H + 
in some places (such as at the left end of the H beam region) and a deficiency in others 
(such as at the conduction front). 

The profile of the electron temperature in the region behind the conduction front 
shows continuous electron cooling throughout the run (Figure 5c). The electrons cool 
through reflection in the conduction front, the energy loss arising from the rightward 
movement of the electric field. This energy loss (electron energy is almost all thermal; bulk 
kinetic energy is a negligible portion of the total) is the source of the free energy driving 
the ambipolar field and, indirectly, the instabilities. The total energy loss of electrons 
behind the shock front is ~ 25% at the end of the run, with the minimum temperature 
being about 50% of the initial electron temperature. 

The bulk velocity and effective temperature of the H + and He + ions is shown in 
Figure 6 for u>j{t = 75, 150, and 300. The effect of the differential acceleration of the two 
ion species is evident in Fig. 6a, in which the heavy ions have a peak bulk velocity which 
is less than half the H + peak bulk velocity. The He + bulk velocity later increases at the 
expense of the H + bulk motion (Fig. 6b) due to the ion-ion instability. Indeed, by the 
end of the run (Fig. 6c), the He + peak kinetic energy is ~ 1.8 times greater than the H Jr 
beam region kinetic energy, for the H + velocity in the beam region has fallen to ~ 0.53 c e 
and the He + peak velocity has risen to ~ 0.36 c e . 

The heating seen in Figures 6d-f is caused by several mechanisms. The H + tempera- 
ture spikes are caused by the presence of fast beam ions, so they are due to an acceleration 
process which creates a secondary beam of ions rather than to a heating process. The 
He + temperature spikes are signatures of the heating caused by ion-ion turbulence. These 
spikes are twice as high as the local H + temperature in Figures 6e and f, and thus a clear 
signature of the preferential heating of the heavy ions. Ion acoustic turbulence is respon- 
sible for the steady rise in H + temperature to the left of the beam region as well as the 
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Fig. 5: Ion density profiles for (a) H + and (b) He + with the density normalized to the initial electron density in Region A. 
(c) Electron temperature profile with the electron temperature normalized to the initial ion temperature. 
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Fig. 6: Ion bulk velocity and temperature profiles. The bulk velocity is normalized to 
the initial value of the acoustic velocity; the temperature to the initial ion 
temperature. 















rise in temperature in the ambient region to the right of the H+ spikes. Both ion acoustic 
and ion-ion turbulence heat the H ions in the beam region and in the fronts. 

If the entire simulation system were to correspond to a soft X-ray source region, 
then distribution functions for the whole system should bear qualitative similarities to soft 
X-ray line profiles. Figures 7a and b shows such distribution functions for (a) H+ and (b) 
Ht+ for the beginning, middle and end of the run. Both species have a positive velocity 
(or blue-shifted) component halfway through the run which grows as a proportion of the 
total distribution function with time, though it never exceeds the stationary component in 
size. The moving component is more substantial; by contrast, the He + component is 
broader and, unlike the H ^ component, increases in velocity with time due to acceleration 
of ions by the ion-ion instability. For neither species does the portion of the stationary 
component for v < 0 broaden significantly during the run. For an electron temperature of 
10 keV, the effective width of the of the He+ distribution (about 200 km/s) and the velocity 
of its blue-shifted component (about 300 km/s), are in rough agreement with velocities 
inferred from soft X-ray observations. Thus, the ion distribution functions reproduce many, 

though not all, of the features seen in soft X-ray line profiles. 

Distribution Functions Across System 



Fig. 7: Ion distribution functions for the entire simulation system. Velocities are normalized to the 
initial thermal velocities of the two species. 

4. CHROMOSPHERIC HEATING 

In this section, the acceleration and heating of ions is examined in a case analogous 
to chromospheric heating. Chromospheric electrons are assumed to be heated to coronal 
temperatures, so the electrom temperature is uniform through the system, as in Section 
3. The electron temperature is elevated with respect to the ion temperature in the chro- 
mosphere, but not in the corona, so vx% is no longer much less than c a throughout the 
system, as is true in Section 3. As a result, a pressure excess in the dense region still 
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drives an outward propagating ambipolar electric field, but no conduction front or density 
wave forms and little ion heating occurs. The density jump of 20 in this section is larger 
than in Section 3, a difference that is not examined in detail because simulations run with 
this density jump but no ion temperature gradient yield essentially the same qualitative 
features as seen in Section 3. The parameters of the simulation are the same as in Section 
3, except n A fn c = 20, T e /Ti A = 100, and T e /T iC = 1. 

The ion phase space diagrams are shown in Figure 8. Figures 8a and b are the 
initial phase space plots of the H + and He + ions respectively. As in the coronal case, 
the source region ions are accelerated outwards by the ambipolar field created by the 
expansion of the heated source region electrons. A difference from the previous simulation 
is that the ambient, or coronal, ions are not swept outward by a conduction front. Instead, 
many of the coronal ions are able to penetrate the ambipolar field region (delimited by the 
dashed lines in Figures 8c-f) because, since the coronal ion thermal velocity Vfic ~ c »t 
they have sufficient kinetic energy to overcome the potential energy barrier the field region 
represents. As a result, no conduction front or shock front forms, and neither does a 
density wave (Figs. 8c, e, and g). 

The region containing the ambipolar field initially confined to the transition region, 
penetrates deeply into both the chromosphere and corona during the run. The left edge of 
the field region propagates to x/A ~ 128 by the end of the run, with an average velocity 
of ~ 0.7 c a . The right edge (indicated by the right side arrows) propagates to x/A ^ 1280 
by wj/i = 300 for an average velocity of ~ 2.5 ( T e /mn J 1 / 2 , which is only slightly greater 
than the local sound speed once the ion temperature is taken into account. While the field 
region expands, the potential drop across it remains roughly constant (~ 4 T e ), so the 
average field strength in the region decreases monotonically with time. 

The initial velocity of many transition region ions matches or exceeds that of the 
outgoing electric field, whose velocity averages ~ 1.2c*. As a result, these ions are acceler- 
ated more efficiently by the ambipolar field, i. e., they gain more energy in the potential 
drop of the field than slower ions, in a fashion similar to the mechanism that produced the 
fast beam ions of the last section. For example, transition region II + ions that begin with 
a velocity of ~ c* in Figure 8a, have been accelerated to a velocity of > 2 c* in Figure 8e, 
an energy gain of > 3 T e . Initially fast He+ ions are likewise accelerated efficiently; ions 
that begin with a velocity of ~ 0.5 c* (Fig. 8b) are accelerated to velocities of ~ 0.8 c a , 
an energy gain of ~ 2 T e . These fast ions form the leading edges of the outfowing ion 
“beams.” 

The effects of the ion-ion streaming instability are visible in the phase space plots of 
both ion species. The velocity dispersion of source region He + ions to the right of ~ x\ is 
increased between u^t — 75 and 150 (Figs. 8d and f) due to the instability, and increases 
further between u>nt — 150 and 300 as ion-ion turbulence grows in strength. This heating 
is accompanied by a weak low energy tail in the H + beam (Fig. 8g). The heating of both 
species is weaker than in the coronal case due to the presence of hot coronal ions which 
effectively decrease the ratio of T e /Ti and thereby suppress the ion-ion instability. 
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Figure 9 shows the bulk velocity and temperature profiles of the chromospheric and 
transition region ions only, since these parameters for the coronal ions remain largely un- 
changed over the course of the run. The ambipolar field more readily accelerates the faster 
transition region ions to high velocitie than the slower transition region and chromospheric 
ions. Such an acceleration hierarchy results in outflowing particle beams with monotoni- 
cally increasing bulk velocity profiles for both ion species (Fig. 9a). The H+ beam is more 
extended than the Hc+ beam because of the greater velocity of the H+ ions. The peak 
kinetic energy of the Hc+ ions initially is also smaller than the H+ peak energy due to the 
more gradual acceleration of He~*~ ions by the ambipolar field (Fig. 9a). As the leading 
edge of the Ht+ ion beam undergoes further acceleration, its velocity becomes half that 
of, and represents the same kinetic energy as the peak ff + velocity (Fig. 9c). 

The heating caused by the ion-ion instability is evident in Figures 9c and d, in which 
the Ht+ temperature exceeds the H+ temperature between x/A ~ 400 and 700 in Figure 
9c, and between x/A ~ 300 and 900 in Figure 9d (the H+ heating actually is due more 
to ion acoustic turbulence than ion-ion turbulence). The extent of the heavy ion heating 
grows between wj/t = 150 and 300, indicating that ion-ion turbulence peaks in strength 
after u^t — 150. 

In order to compare these results to soft X-ray line emissions, we examine the 
distribution functions of the two ion species, which are shown in Figure 10 (tf + in Fig. 
10a, Ht+ in Fig. 10b). As seen previously in the coronal case (Fig. 7), the blue-shifted 
component of both ion species increases in size as the run continues, with that of H+ 
consisting of a greater proportion of the total distribution function than the blue-shifted 
component of Ht+ at any time t > 0. The central difference is that the velocity of the peak 
of the blue-shifted component increase with time in both distributions rather than just the 
Ht+ distribution due to the continuous acceleration over time of source region ions by 
the ambipolar field. The effective temperature, or velocity dispersion, of the blue shifted 
components of both species also appears to increase with time, to a value approaching T e 
by the end of the run. The broadened distributions bear close qualitative similarities to 
the distribution functions for the coronal heating simulation (Fig. 7). However, if one 
assumes that the electron temperature is 100 eV, the effective thermal velocity of the 
distribution, as well as the velocity of the moving component, is only about 50 km/s, which 
is too low to explain the observations, however. 

5. CONCLUSION 

Ion heating associated with the expansion of hot electron populations in both coro- 
nal chromospheric applications during solar flares have been analyzed using particle simu- 
lations as methods of reproducing heavy ion soft X-ray line features observed during solar 
flares; specifically, nonthermal line broadening and the presence of blue-shifted compo- 
nents. In the simulations, the expansion of heated electrons across density gradients is 
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Fig. 10: Distribution functions for all the ions in the simulation system. Velocities 
are normalized to initial chromospheric thermal velocities. 
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found to create an ambipolar electric field which is capable of accelerating ions up to ener- 
gies of order the lectron thermal energy. In the coronal case, the H ions form a constant 
velocity and density beam region between forward and backward propagating fronts of the 
ambipolar field, the forward front eventually steepening into an electrostatic shock. The 
ions in the beam region are unstable to the ion-ion streaming instability, which then acts to 
preferentially heat and accelerate the jions by extracting the H+ beam energy. This 
nonthermal heating as well as the acceleration is clear in the ion distribution functions, 
which bear some similarity to soft X-ray line profiles* In the chromospheric case, no similar 
beam region is formed, the ion-ion instability is weaker, and little preferential heavy ion 
heating occurs. However, substantial broadening is present in the ion distribution func- 
tions due to the nonuniform motion of the ions. Thus, these simulations have been able to 
reproduce some of the features, though not all, of soft X-ray line emissions. Only in the 
coronal heating simulation are the broadening and the speed of the blue-shifted component 
in agreement with soft X-ray line observations. 

Future work on this project includes the addition of a heavy ion with a charge 
to mass ratio of about 2 to represent those ions actually emitting the soft X-rays, and 
the addition of a method to simulate the continuous heating of the system by energetic 
electrons. 
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